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ABSTRACT 

It has been suggested that the 170 day period in the light curve of the low mass 
X-ray binary 4U 1820-30 arises from the presence of a third body with a large incli- 
nation to the binary orbit. We show that this long period motion arises if the system is 
librating around the stable fixed point in a Kozai resonance. We demonstrate that mass 
transfer drives the system toward this fixed point, and calculate, both analytically and 
via numerical integrations, that the period of libration is of order 170 days when the 
mutual inclination is near the Kozai critical value. The non-zero eccentricity of the 
binary, combined with tidal dissipation, implies that the rate of change of the binary 
period would be slower than, or even of opposite sign to, that implied by standard 
mass transfer models. If the 170 day period results from libration, then, contrary to 
appearances, the orbital period of the inner binary is increasing with time; in that case, 
{e/0.009fQ/k2 > 2.5 x 10^ where k2 ^ 0.01 is the tidal Love number and e = 0.009 is 
the fiducial eccentricity of the inner binary. It appears unlikely that the observed neg- 
ative period derivative results from the smaller than expected (but positive) value of P 
combined with the previously suggested acceleration of the system in the gravitational 
field of the host globular cluster NGC 6624. The discrepancy between the observed 
and expected period derivative requires further investigation. 

Subject headings: binaries: LMXB — stars: individual 4U 1820-30 — stellar dynamics — 
celestial mechanics 

1. INTRODUCTION 

4U 1820-30 is a low mass X-ray binary (LMXB) located near the center of the globular cluster 
NGC 6624. The binary orbital period is Pi ~ 685 s, revealed in X-ray observations as a modulation 
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with ~ 2-3% peak to peak amplitude (|Stella et alJll987|) . Subsequently. lAnderson et alj (Il997|) 
discovered a ~ 16% peak to peak modulation (period 687.6 ± 2.4 s) in the UV band from HST. 

This short period, lo w amplitude variation is very stable, with /'//' = (-3.47 ± 1 .48) x I0~^yr~^ 
Jchou & GrindlavlboOll). which i s consi stent with the earli er measurement of P/P = (-5.3 ± 1 . 1) x 
I0~^yr~^ from Ivan der Klis et alJ (|l993ar) : this stability led lChou & Grindlayl (|200lh to suggestion 
that this modulation reflects the orbital period of the binary. 

Both the short binary period and the type I X-ray bursts observed in this system imply that 
the secondary star is a helium white dwarf, of mass m2 = (0.05-0.08)A/q, accreting mass onto a 
primary neutron star (IRappaport et al.ll 19871) . The distance to the source is estimated to be 7.6 ± 



0.4kpc (IKuulkers et al 



20031 



It is striking that neither the magnitude nor the sign of the period derivative is consistent 
with the prediction P/P > x 10~^yr~' of the standard evolution scenario for compact binaries 
overflowing their Roche lobe feappaport et al. 1987). It has been suggested that the negative 
period derivative is only apparent, i.e., that it is not intrinsic to the binary, but instead reflects 
the acceler ation of the binary in the gravitational potential of the globular cluster which houses 
the binary (Ivan der Klis et al.lll993bl) . However, quantitative estimates show that the acceleration, 
while of roughly the right magnitude, is unlikely to be large enough, by itself, to explain the 
large discrepancy between the evolution scenario and the observations dvan der Klis et al.lll993bl : 



King etal.ll 19931 : IChou & Grindlavll200ll) . 



A second striking property of 4U 1820-303 is the much larger luminosity variation, by factor 
of > 2, seen at a period of /3 ~ 171 days. Analysis of the RXTE ASM data shows that this long pe- 
riod m odulation does not exhibit a significant period derivative, P3/P3 < 2.2 x lO'^^vr'^ llChou & Grindlay 
200 ll) . The ratio between this long period and the binary orbital period is ~ 2 x 10"^, which ap- 



pears to be too high to be due to disk precession at the mass ratio of the system (lLarwoodlll998 



Wijers&Pringlell 19991) . 



In this paper we adopt the assumption of iGrindlayl (Il988l) . that the 171 day period is due to 
the presence of a third body in the system. The third (outer) star modulates the eccentricity of the 
binary at long term period P3 ~ /'|/(e/'i), where Pj is the orbital period of the third star and e is 
the eccentricity of the inner binary. Taking into account only perturbations from the third star, the 
binary orbital period of 685 s and ~ 171 day long-term modulation imply that the orbital period 
of the third star must be ~ 1 day. The presence of additional sources of precession, such as that 
due to tidal distortion of the white dwarf secondary, requires a stronger perturbation from the third 
body and hence a smaller orbit in order to modulate eccentricity of the inner binary at the 171 day 
period. We show that the luminosity modulation arises from variations in the eccentricity of the 
inner binary associated with libration around a stable fixed point in the Kozai resonance. 

Tidal dissipation in the white dwarf, driven by the eccentricity of the binary orbit, tends to de- 
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crease both the eccentricity and the semimajor axis (hence period) of the binary, which we suggest 



is res ponsible, in part, for the anomalous observed period derivative — note that Rappa port et al 



(fl987 ) did not treat the effects of tidal dissipation. The combination of tidal dissipation and mass 
transfer will result in a lower value of P/P than that produced by conservative mass transfer alone. 

For rapid enough dissipation, or, expressed another way, for low enough values of the tidal 
dissipation parameter Q, P <0 could result. We do not favor this as the explanation for the ob- 
served negative period derivative; we show that such rapid dissipation damps eccentricity within 
lO""* of the system's lifetime. Subsequently the mass transfer takes over the evolution of the semi- 
major axis. In other words, we would be incredibly lucky to observe the system in the short time 
that e is significant, in the absence of another perturbing influence. We also show that, given 
the most recent estimates for the acceleration of millisecond pulsars in the gravitational field of 
the globular cluster, the cluster gravity does not appear to contribute significantly to the observed 
period derivative of 1820-30. 

Thus it appears that, while both tidal dissipation and acceleration in the gravitational field of 
the cluster contribute negatively to the period derivative, they can not fully explain it. Since we 
favor the hierarchical triple model as an explanation for the origin of IVlday period of luminosity 
variations, we suggest that the apparent negative period derivative, which is a 2 - a result, may 
either be an observational artifact or due to the some yet not understood physical processes. 

The relation between the luminosity variations and the period derivative is deeper; we argue 
that the (intrinsic) increase in the semimajor axis of the binary (driven by Roche lobe overflow) 
leads to trapping of the system deep in the Kozai resonance. The resonance transfers angular 
momentum from the inner binary to the third star, and back, periodically, without affecting the 
semimajor axis of either orbit. However, the dissipation associated with the strong tides when the 
forced eccentricity is largest does remove energy from the orbit of the inner binary. This energy loss 
peaks when the mutual inclination is small. It is well known that this coupled Kozai-tidal evolution 
tends to leave the system with a mutual inc lination between the two orbits near the Kozai critical 
value (~ 40°); see, for example. Figure 4 in IWu et al.l (|2007r) or Figure 7 in lFabrycky & Tremaine 
(|2007|) . We show that the period of small oscillations is naturally ~ 170 days when the mutual 
inclination is close to the Kozai critical value. Whether the evolution of the inclination in systems 
like 1820-30, which, unlike the planetary systems, is know to undergo Roche lobe overflow, is a 
question we are currently investigating. 

This paper is organized as follows. In Owe develop an analytic understanding of the system, 
describing the resonance dynamics, calculating the location of the fixed point as a function of the 
system parameters (stellar masses, orbital radii, and the mutual inclination of the two orbits), and 
the frequency (or period) of small oscillations. In Owe describe a possible dynamical path by 
which the system arrived at its present configuration. The dynamical history relies crucially on 
both the Roche lobe overflow (which drives the system into resonance) and the tidal dissipation. 
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which tends to drive the mutual inclination toward the Kozai critical value. In section|4]we describe 
the results of numerical integrations of the equations of motion, presenting a fiducial model that 
reproduces the observed properties of 4U 1820-30. We also demonstrate trapping in the case of 
an expanding inner binary orbit, and detrapping in the case of a shrinking binary orbit. In ^ we 
use the model to put constraints on the ratio of the tidal dissipation parameter Q and the tidal Love 
number (^2) of the Helium white dwarf for our fiducial eccentricity. We discuss our results, and 
those of previous workers, in ^ We present our conclusion in the final section. We give the details 
of the numerical model in the appendix A. In appendix B we discuss in details adiabatic invariance 
of the action and how it governs the evolution of the system by comparing analytic and numerical 
analysis. 



2. UNDERSTANDING THE DYNAMICS OF THE 4t/ 1 820 - 30 SYSTEM 

The presence of a third body orbiting the center of mass of a tight binary will induce changes 
in the orbital elements of the binary, changes that take place over a variety of tim e scales. The 



changes are particularly dramatic if the mutual inclination of the two orbits is large. iKozail (I1962|) 
showed that when the initial inclination between inner and outer orbits has values between some 
critical inclination ic-it and I'&O" -lent, both the eccentricity of the inner binary and the mutual 
inclination undergo periodic oscillations known as Kozai cycles. 

The period of the Kozai cycles is much longer than either the binary's orbital period, or the 
period of the outer orbit. This justifies the use of the secular approximation, which involves aver- 
aging the equations of motion over the orbital periods of inner and outer binaries; as a result, the 
averaged equations of motion predict that the semimajor axes of both binaries are unchanged. 

If the luminosity variations in 4U 1820-30 are due to Kozai cycles, the semimajor-axis ratio 
ii-out/o. ~ 8, so in our analytic work we use the quadrupole approximation for the potential experi- 
enced by the inner binary due to the third body. In our numerical work we keep terms to octupole 
order, but we show that the higher order terms change the quantitative results only slightly. 

The angular momentum of the outer binary is much greater than that of the inner, so that the 
orientation of the outer binary is, to a good approximation, also a constant of the motion. In that 
case, after the averaging procedure, the final Hamiltonian has one degree of freedom. 

Kozai cycles are the consequence of a 1 : 1 resonance between the precession rates of the 
longitude of the ascending node Vt and the longitude of the periastron (D of the inner binary. The 
condition for Kozai resonance, cij-fi = 0, is satisfied only for high inclination orbits; for low incli- 
nations, the line of nodes precesses in a retrograde sense {(l < 0), while the apsidal line precesses 
in a prograde sense. 

We employ Delaunay variables to describe the motion of the inner binary. The angular vari- 
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ables are the mean anomaly I, the argument of periastron lu, and the longitude of the ascending node 
fi; of these, only u appears in the averaged Hamiltonian. Their respective conjugate momenta are: 



C 

Q 
U 



Ga 



m\m2 

= CVl-e^ 
= Qcosi. 



(1) 

(2) 
(3) 



The longitude of periastron is (D = + Recall that we are assuming that the semimajor axis of 
the outer binary is large enough that the total angular momentum is dominated by that of the outer 
binary, so that i is effectively the mutual inclination between the two binary orbits. We occasionally 
refer to the elements of the third star, using a subscript 'out' to distinguish them from those of the 
inner binary. 

After averaging over / and lout, the Hamiltonian describing the motion of a tight binary orbited 
by a third body, allowing for the effects of both tidal and rotation al bulges on the secondary, and 



for the apsidal precession induced by general relativistic effects, is dlnnanen et al.ll 19971 : iFord et al. 
200d : lFabrvckv & Tremain3l2007h 
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where the term proportional to A is the Kozai term, the term proportional to B enforces the average 
apsidal precession due to general relativity, and the terms proportional to C and D represent the 
tidal and rotational bulges, respectively; the explicit appearance of the tidal Love number ^2 in the 
latter two terms highlights the fact that these terms represent the effects of the white dwarf's tidal 
and rotational bulges. The expressions for the constants are 
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Recall that the semimajor axis and eccentricity of the outer body's orbit are denoted by aout and 
eouf The quantity = 2Gm\/c^ in equation Q is the Schwarzschild radius of the neutron star. 

As just noted, the term proportional to D accounts for the rotational bulge produced by the 
spin of the white dwarf. The spin is projected onto the triad defined by the Laplace-Runge-Lenz 
vector, pointing along the apsidal line from the white dwarf at apoapse toward the neutron star, and 
denoted by a subscript e, the total angular momentum vector, subscript h, and their cross product, 
denoted by q. We have scaled the spin to the orbital frequency (or mean motion) n, so that, e.g., 
(le = VLe/n. We do so because we anticipate that for small eccentricity the white dwarf will be 
tidally locked. Then /((^spin) = '2^h~^l~^l ^ dimensionless quantity of order unity. 

For the fiducial values of the system parameters listed in table 1, A ^ 1.73 x 10'*'*, the ratios 
B/A ^ 0.53, C/A ^ 1.82, and D/A ^ 2.54. 



2.1. The Kozai mechanism 



We start our discussion of the dynamics of the system by focusing on understanding the Kozai 
mechanism, neglecting forces due to the tidal and rotational bulges of the Helium white dwarf in 
the inner binary, and the effects of general relativity. 

We locate the resonance by looking for a fixed point of the Hamiltonian; since we are neglect- 
ing the tidal and rotational bulges, and the general relativistic precession, we set 5 = C = D = and 
differentiate the Hamiltonian with respect to uj, to find ujf = 0, 90°, 180°, 270°. The fixed points at 
Uf = 90° and ujf = 270° are stable. Differentiating the Hamiltonian with respect to Q, substituting 
CO = 90° (or 270°) and setting the result equal to zero, we find Qj = {5/3)'H^ C^. In terms of the 
eccentricity, 

ef= ^~[-^coshf, (10) 

where the subscript / indicates that this is the eccentricity of the stable fixed point. The frequency 
of small oscillations around the fixed point (small librations) is 



■'f'Sf 



1/2 



(11) 



Performing the derivatives. 



uo = uja \ 18 + 90- 




where we have defined 
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From equation (flOl) we see that the critical inclination for a Kozai resonance to occur, in the 
absence of other dynamical effects, is i^it = cos"^ a/3/5 ^ 39.2°. If i > idt, orbits started at w = 90° 
with e < Cf will librate around the fixed point, so that u remains between 0° and 180° (or an even 
more restricted range). From equation (fT2l) or (fT4l) . the period of small oscillations Pq ~ 1/e/, a 
point that will be important later. 

In contrast, orbits started at w = and e > will circulate {u will range from to 360°). 
Librating and circulating orbits are separated by the separatrix, an orbit that neither librates nor 
circulates. The width of the separatrix (as measured by the excursion in e) depends only on the 
initial inclination: e^ep = [1 - (5/3) cos^ z]^/^. 

Examples of librating and circulating orbits (for a system including the effects of GR and tidal 
bulges) are shown in section |4l 

Note that, even for systems with i < i^it, where no stable Kozai fixed point exists, both the 
mutual inclination and the eccentricity of the inner binary can undergo oscillations with significant 
amplitude (although reduced compared to the case with i > lent)- 

Kozai cycles will be substantial only as long as the perturbation from the outer body dominates 
over the other sources of apsidal precession in the inner binary orbit, a point we now address. 



The physical effects represented by the terms proportional to B, C, and D are capable of 
suppressing Kozai oscillations. We investigate their effects in this section. 

As an aside, there is a small apsidal precession introduced by dissipative effects in the He 
white dwarf, but this precession rate is negligible compared to the other three. We mention it here 
because tidal dissipation has a major role to play in the capture (or otherwise) of the system into 
the Kozai resonance. 

The equations for the precession rates due to the Kozai mechanism, that due to general rela- 
tivity, and the tidal and rotational bulges of the white dwarf, are: 



2.2. Kozai cycles in the presence of additional forces 
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The Kozai term (equation [T5] ) can be either positive or negative, depending on the value of 
sinj. Both the white dwarf tidal bulge and the GR terms are positive, so both tend to suppress 
Kozai oscillations. The term induced by the white dwarf rotational bulge, on the other hand, can 
be of either sign, depending on the orientation of the white dwarf spin. If the white dwarf is tidally 
locked and if its spin is aligned (which we assume in our analytic model, but not in our numerical 
models), this term contributes positive u. In case of non-aligned spins the precession rate may be 
negative (as we will see). 



2.2.1. The tidal bulge and the tidal Love number ^2 

The tidal bulge of the white dwarf in 4U 1820-30 dominates the non-Kozai apsidal precession 
rate, for physically plausible values of ^2- We remind the reader that in Newtonian gravitational 
theory the tidal Love number ^2 is a dimensionless constant that relates the mass multipole moment 
created by tidal forces on a spherical celestial body to the gravitational tidal field in which it is 
immersed; in other words, ^2 encodes information about body's internal structureQ 

We use k2 = 0.01, which is computed by Arras (private communication) as the ratio of the 
potential due to the perturbed mass distribution, to the external potential causing the perturbed 
mass, under the assumption that our He white dwarf is a fluid object. 

Soft X-ray observations of the source indicate a rather small absorption, consistent with that 
expected to be produced by the interstellar medium of the Galaxy; this rules out any significant 
outflows from the accretion disk or the surface of the white dwarf. This implies an absence of 
mass loss through the L2 Lagrangian point o f the white dwarf, w hich puts an upper limit on the 



eccentricity of the inner binary; according to iRegos et al.l (|2005|) . for our system parameters, the 



upper limit on the eccentricity of inner binary is e^ax — 0.07. 



'in a confusing usage, the apsidal precession constant, which is a factor of two smaller than the tidal Love number, 
but which we do not employ, is also denoted by ^2- 



9 



If 4U 1820-30 has a non-zero but small eccentricity, as indicated by the observed luminosity 
variations, then in the absence of a third body, the precession rate of the binary orbit is dominated 
by the tidal bulge induced in the white dwarf by the gravity of the neutron star; from equations 
(fT6l ) and (flTI ), the tidal bulge induces a precession rate at least a few times that induced by GR: 



In order for the Kozai mechanism to produce significant variations in e, the Kozai-induced 
precession rate must be comparable to or larger than the sum of the precession rates produced by 
the other terms. For physically realistic values of k2, as we have just seen, the precession rate 
induced by the tidal bulge of the white dwarf is by far the largest, so if the Kozai effect is to be 
important, it must produce a precession rate larger than ujjb- 

2.3. Libration around the fixed point and the frequency of small oscillations 

2.3.1. Why libration? 

For the values of the tidal Love number and eccentricity listed in table 1, the period of the 
precession rate induced by the tidal bulge, Pjb = 27r/a)TB, is a factor of ten shorter that the period 
of the observed luminosity variations. If this term set the rate of precession, and the eccentricity 
varied as a result of this precession, then the variations in X-ray luminosity would occur with a 
period substantially shorter than the observed 170 days. 

In order to produce a much longer period, some other term must tend to produce a negative 
precession rate. When this negative precession rate is added to that produced by the tidal bulge, 
the resulting period can be much longer than that produced by the tidal bulge alone. 

Under the assumption that the white dwarf is tidally locked (we show later it is not), the only 
term capable of producing a negative precession rate is the Kozai term. Hence we are led to look 
for a cancellation between the Kozai precession rate and the precession rate induced by the tidal 
bulge. 

However, it is not enough to ask for a rough cancellation. To get the observed precession rate, 
the sum of all the terms must cancel to better than 10%. This requires some fine tuning of the 
mutual inclination, a rather unsatisfactory situation. 

On the other hand, if the system is captured into libration, then the sum of all the precession 
terms is exactly zero. If the system is deep in the resonance, then the period of libration is simply 
the period associated with small oscillations around the fixed point. We show here that the period 
of small oscillations is naturally around 170 days, if the mutual inclination is near the critical value 
for Kozai oscillations. 




(19) 
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2.3.2. The frequency of small oscillations 

Setting the first derivative of the Hamiltonian (Hj) with respect to tu and Q to zero, we find the 
following expression for the location of the stable fixed points in the limit of small eccentricity: 

Uf = 90° , 270° (20) 

(21) 



^/ = 

We can write the second of these as 



18-30^-f-120yt2f-3)t2f 
\\ 60f + |f + 8m2f + f^.f ■ 



30 cos^Zm, -cos^z 



60f + |f + 840fc2f + f^2f' 



(22) 



where 



2 3 1 5 C 1 D 

cos z„,Y = -Ako — - — kj — . (23) 

5 30A A 10 A ^ 



Evaluating the second derivative of the Hamiltonian at the fixed point we obtain the expression 
for the frequency of small oscillation around the fixed point: 



-1/2 r2\ R r3 / rll r9 rl 

18+90^^ +2-±- + fc2- 3150^-1680^ + 90^ 
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1/2 



g/sinz/, (24) 



which should be compared to equation (fT2l) . As in the pure Kozai case, the period of small oscil- 
lations Pq ~ 1/^/- 

Figure [U shows Pq as a function of the initial inclination. As the initial inclination increases 
above the critical value, the period of small oscillations decreases rapidly. Increasing the initial 
inclination increases the magnitude of the Kozai torque; in the absence of other torques, and for 
inclinations above the critical inclination, increasing the magnitude of the Kozai torque is analo- 
gous to increasing the restoring force in a harmonic oscillator, thereby increasing the frequency of 
oscillation. When there are other torques in the problem, the critical inclination will change; for 
example, the presence of a tidal bulge on the secondary increases the critical inclination. 

Very near the critical inclination, the effective restoring force is small, ~ sinz/, so the 
frequency of small oscillations is small, and the period of oscillations is large — hence the rapid 
increase in Pq as the inclination decreases toward the critical inclination {i^rit ~ 44° in Figure [T]). 

Figure [2] shows /q as a function of aout- As expected from the nout ~ a^ut dependance of ua, 
the period of eccentricity oscillations increases rather rapidly with aout- 
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Fig. 1. — The period of small oscillations vs. the initial inclination for a system similar to 4U 
1820-30, with ^2 = 0.01. The critical inclination is icru ~ 44.7°. At large inclinations, well above 
icrit, the period of libration is of order days. Only if i ^ z,„y is the period of order 170 days. The 
solid line is the prediction of equation (l24l) : the solid circles come from numerical integration of the 
equations of motion to quadrupole order, while the open squares come from integration accurate 
to octupole order. 
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Fig. 2. — The period of small oscillations vs. Uout- As the semimajor axis of the outer binary, Qout 
increases, the period of small oscillations is increasing too, which is expected from the Uout ~ <^out 
dependance of coKozai- The solid line is the prediction of equation (|24|) . while the solid circles and 
open squares are from numerical integrations accurate to quadrupole and octupole order respec- 
tively. 
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3. MASS TRANSFER, TIDAL DISSIPATION, AND CAPTURE INTO LIBRATION 

We have shown that physically plausible values of ^2 lead to a precession frequency ujtb that is 
much larger than the observed frequency of luminosity variations in 4U 1820-30. We then appealed 
to an equally large precession, of the opposite sign, produced by the Kozai interaction, to cancel 
the prograde precession caused by the tidal bulge. In order to avoid fine tuning, we argued that 
the system has to be in libration, so that the observed low frequency actually arises from libration, 
rather than precession of the apsidal line of the binary orbit. 

Whether the tidal bulge or GR effects produce a larger precession rate, we argue that it is 
no coincidence that the magnitude of the Kozai precession rate is equal to the sum of the other 
precession rates: the system will evolve so as to capture the orbit into resonance, in which the sum 
of all the precession rates is zero. 

Capture into libration in the Kozai resonance is a natural consequence of semimajor axis 
expansion, the latter driven by mass loss from the white dwarf as a result of its overflowing its 
Roche lobe. The action J Qdu is an adiabatic invariant (for detailed discussion see appendix B), 
since the semimajor axis of the binary orbit is expanding on the accretion time scale m2/m2 ~ 
10^ yr, much greater than either the orbital or precession time scale. In contrast to mass transfer, 
tidal dissipation tends to shrink the semi-major axis; if this effect dominates, trapping into the 
Kozai resonance is not possible. 

How does expansion of the inner orbit lead to capture into libration? As a increases, the mu- 
tual torque between the two orbits will increase as well — the inner orbit is expanding, effectively 
moving closer to the outer orbit. This increasing torque corresponds to a deepening of the Kozai 
potential, and an expansion in the size of the separatrix of the Kozai resonance. Orbits other than 
the separatrix have a fixed action, while the action of the separatrix is increasing. If the increase 
in the action of the separatrix grows to exceed the action of an initially circulating orbit, that cir- 
culating orbit will be captured into resonance, and begin to librate. As a continues to expand, the 
captured orbit will move closer and closer to the fixed point of the resonance, librating with the 
frequency of small oscillations. 



More quantitatively, mass transfer tends to increase a (IRappaport et al.lll982() : 
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0.4242 



a 



9/2 



where q = mi/m2 and K = kOy / (unip); k is the Boltzmann constant, /i is the mean molecular weight, 
nip is the mass of the proton and 9y is the polytropic temperature. The parameter K is given by the 
following mass -radius relation: 

K = N,Gml-''^"^R^'^"^-\ (26) 



where A^^ is a tabulated numerical coefficient (for n = 1.5 it is 0.4242: IChandrasekhan (|1939() '). 
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Tidal dissipation in the white dwarf will tend to reduce the semimajor axis of the binary. In 
the limit of small eccentricity, 



We argue that the orbit must be expanding, (e/e)™ is 100 times shorter than {ajd)jr,, so 
unless something excites e (such as third body or thermal tides) we are unlikely to catch the system 
in a phase where periastron, = a(l -e), is increasing while a is decreasing. 



Our numerical model treats the gravitational effects of the third body in the quadrupole ap- 
proximation. We average over the orbital periods of both the inner binary and the outer companion. 
We demonstrate in section 14.41 and in figures [T] and [2] that treating the effects of the third body in 
octupole approximation does not qualitatively change our findings. We include the following dy- 
namical effects: 

• periastron advance due to general relativity; 

• periastron advance arising from quadrupole distortions of the helium white dwarf due to both 
tides and rotation; 

• orbital decay due to tidal dissipation in the white dwarf; 

• loss of binary orbital angular momentum due to gravitational radiation; 

• conservative mass transfer from the helium white dwarf to the neutron star primary driven 
by the emission of gravitational radiation. 

Note that the Kozai mechanism described in the previous section is included in the three body 
gravitational dynamics. The equations used in our model are listed in the appendix. 



We use as fiducial parameters mi = 1.4M0, m2 = 0.067Mq, and = 0.55Mq. The semimajor 
axis of the inner binary is a = 1 .32 x 10^° cm, chosen to match the observed orbital period of 685 s. 




(27) 



4. NUMERICAL RESULTS 



4.1. Numerical model using the quadrupole approximation 



4.2. Results 
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The radius of the Helium white dwarf is R2 = 2.2 x 10^ cm, while the fiducial Love number is 

k2 = om. 

To reproduce the 171 day eccentricity oscillations (Figure [3]), we use the following initial 
parameters: Uout = 8.0a (yielding Pout = 0. 15 day). We start with = 0.009, = 90°, = 44.715°, 
and eoutfl = 10""*. 

Figure [3] shows the eccentricity oscillations of the inner binary, with a period of 17 1 days, over 
a decade. The amplitude of the eccentricity oscillations is of order of 7 x 10"^, which is sufficient 
to enhance mass trans fer enough to produce the observed luminosity oscillations of a factor of > 2 
dZdziarski et al.ll2007|) : see their Figure 3. The amplitude of the eccentricity oscillations depends 
on the initial eccentricity, as illustrated in Figured a lower initial eccentricity produces eccentric- 
ity oscillations with higher amplitude. If the system circulates, the amplitude of the eccentricity 
oscillations is larger still. 

Having the system trapped in libration about the fixed point explains both the origin of the 171 
day period luminosity variations, as well as the small amplitude of the eccentricity oscillations; the 
observations require that magnitude of the eccentricity oscillations be small so as to avoid overly 
large luminosity variations — a point we return to below. 



4.3. Resonant trapping and detrapping of 4U 1820-30 



The mass transfer rate is determined by the inner binary mass and semimaior axis. These 



Anderson et al. 



parameters are reasonably w ell constrained from observations (iStella et al.ll 19871 : 
19971 : iRappaport et al.lll987b . The amount of tidal dissipation is parameterized by the tidal dis- 
sipation factor Q, which for white dwarfs is not well constrained at all. If we know the value of 
period derivative, P, we can constrain Q (or more precisely, (e/0.009)^Q/k2, see equation [27]) for 
the white dwarf in the system. 

We argued at the end of O that the intrinsic P must be positive, since a shrinking binary 
orbit and a decaying eccentricity quickly lead to mass transfer driving expansion of the binary 
orbit. There is a second argument against an intrinsic negative P: if the orbit of the inner binary is 
shrinking, an initially librating orbit will quickly become circulating, and the period of luminosity 
variations will change dramatically. If we tune Q to the value that reproduces the observed negative 
period derivative (Q = 2.5 x 10^, assuming ^2 = 0.01) and let the system evolve, the system is driven 
out of libration after about 1500yr, as shown in Figure[5l As the figure shows, the eccentricity of 
the inner binary decreases significantly due to tidal dissipation, which in turn reduces the strength 
of tidal dissipation. With tidal dissipation weakened, mass transfer will dominate the evolution of 
the semimajor axis and, as expected from the standard evolutionary scenario, the semimajor axis 
starts to expand (not shown in the figure). As long as there is some small eccentricity in the inner 
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Fig. 3. — The eccentricity as a function of time (upper panel) and the phase space (e versus u) for 
our fiducial model. The period of the eccentricity oscillations is 171 days, and the amplitude of the 
eccentricity oscillation is sufficient to produce the observed factor of 2-3 variation in luminosity. 
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Fig. 4. — Phase portrait for four different initial eccentricities at initial inclination i = 44.715" and 
initial tu = 90°. Our fiducial orbit, with a libration period of 171 days and the amplitude of the 
eccentricity oscillations sufficient to produce observed variations in the luminosity, is labeled as 
the "current orbit in libration"; it is the same orbit presented on Figure [3l The unlabelled librating 
orbit is very near the fixed point; it has a period of the eccentricity oscillations that is close to 
but shorter than the observed period, while the variations induced in the luminosity are too small 
compared to the observations. The circulating orbit produces luminosity variations that are too 
large, as well as having an oscillation period that is too long. 
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binary there is some tidal dissipation present that tends to slow down the expansion rate of the 
semimajor axis. 

The reason for the detrapping is rather subtle. First, we note that the decrease in e is not due 
to direct tidal damping; Equation (|38l ) predicts (e/e)™ ~ 10^ yr, while e changes by factor of 2 in 
2000 yr. To verify this, we have set eVo = 0^ and verified that integration yields the same result. 
The reason for such a short time scale for decrease in e lies in the fact that the spins do not remain 
tidally locked throughout the evolution of the system and the evolution of the eccentricity is rather 
strongly influenced by the their lack of pseudo-synchronism. Detailed discussion and figures are 
given in appendix B. 

On the other hand, if the observed negative period is not an intrinsic property of the system, 
in other words, if the effect of mass transfer wins over the effect of tidal dissipation, the action of 
the separatrix increases with time, and trapping will occur. 

Figure [6] shows a system initially put on a circulating orbit. As the integration proceeds, 
the separatrix expands, eventually capturing the orbit, which then librates for the duration of the 
integration. 

4.4. Numerical model using octupole approximation 

In this subsection we treat gravitational effects of the third body in the octupole approxima- 
tion. As in the case of the quadrupole approximation, we derive our equations of motion from 
the double averaged Hamiltonian and we include all of the previously listed dynamical effects. As 
Figure|7]demonstrates, the octupole approximation does not change qualitatively our previous find- 
ings. All parameters, except the initial inclination, used in the octupole approximation are listed in 
table 1. In order to produce the 171 days period of the eccentricity oscillations, and the amplitude 
of the eccentricity oscillation that produces the observed factor of 2 - 3 variation in luminosity, a 
slightly higher inclination is required {i = 45.1"). 

5. ON THE VALUE OF Q AND THE ORIGIN OF THE SMALL (OR NEGATIVE) P 

The standard theory of Roche lobe overflow predicts P/P > x 10"^ yr"^ The measured 
P/P = (-3.47 ± 1.48) X 10"*^yr"^ is eight standard deviations away from this value. We have argued 
in previous section that P/P should be positive, but even if it is two or three standard deviation from 
the measured value, it is still five below the predicted value. The origin of this discrepancy has been 
a puzzle since it was discovered. 

The suggestion that the binary has a finite eccentricity immediately suggests a reason for the 
low value of P: tidal dissipation in the white dwarf will tend to reduce the semimajor axis of the 
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Fig. 5. — a) (jj vs t. We start the evolution of the system by placing the system in libration. 
Here we use Q = 2.5 x 10^, which is the value required to reproduce the observed negative period 
derivative. The action of the separatrix is decreasing, and the system is ejected from the resonance 
after about 1500yr. b) e vs to, phase space evolution plot, showing that the orbit evolves from 
libration to rotation, with the transition occurring between the 1500 and IVOOyr snapshots, c) P/P 
vs t. The period derivative remains negative only for about IVOOyr, which is 10"^ of the lifetime 
of the system; the eccentricity damps sufficiently that the m term takes over. 
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libration. Here we have used 2 = 5 x 10^, while the initial eccentricity is eo = 0.032; all other 
parameters are the same as used in Figured b) e vs a; for the same integration. 
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Fig. 7. — The eccentricity as a function of time (upper panel) and the phase space (e versus a;) in 
the octupole approximation. To produce the 171 days period of the eccentricity oscillations, and 
the amplitude of the eccentricity oscillation that is sufficient to produce the observed factor of 2-3 
variation in luminosity, slightly higher inclination is required (« = 45. 1"). 
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orbit, contributing a substantial negative term to P. 

The tidal dissipation could in fact dominate the orbital evolution, overcoming the effects of 
mass transfer as seen in Figure[51 We do not argue for this point of view, however, because it would 
be unlikely that the system could be observed in a stage of the evolution that last only 10"^ of its 
lifetime. In addition, we believe that the system is trapped in libration. 

The observed P/P consists of at least three parts: 

L) jt) JL) J^:) . (28) 

^ J ohs \P J Roche J accel J TD 

The values of the observed and Roche terms were given above, and, as noted there, they are not 
consistent with each other. The second term on the right hand side of equation (|28l) represents the 
acceleration experienced by 1820-30 in the gravitational field of its host globular cluster, while the 
third term on the right represents the effects of tidal dissipation in the white dwarf secondary. 

A natural explanation for the observed negative P might be provided by a combination of 
the last two effects, but still allow for the system to be trapped in resonance. First, tidal dissipa- 
tion reduces the intrinsic P/P substantially from that expecte d due to Roche lobe overfl ow alone. 



but leaves P/P > 0. We then we appeal to the argument of Ivan der Klis et al.l (Il993b|) . that the 



(P/P)accei term produces an apparent ne gative total P. Indeed, g iven the most recent published 



estimate of a,„ax/c = 7.9 x 10 ^yr ' from van der Klis et al.l (| 1993b) . it is plausible that we would 



observe a negative period derivative, while the intrinsic (or physical) period derivative is in fact 
postive. 

However, recent estimates for the cluster acceleration from millisecond pulsar timing suggest 
a maximum of a,nax/c = 1.3 x 10"^jr"' ( Lynch and Ransorn , private communication), an order of 



magnitude smaller than the estimate from lvan der Klis et al.l (Il993b|) : if the smaller value holds up. 



the observed negative period derivative is difficult to understand in the context of current models. 

Given that the measured negative period dervative is significant only at the two-sigma level, 
and that there is no clear physical explanation for such an orbital decay, it is worth consider- 
ing the possibility that the observed value is in error. If we ignore the observed negative period 
derivative, and simply assume that the intrinsic P is postive, we find a lower limit on Q given by 
(e/0. 009)^2/^2 > 3.15 X 10^. We can get a firmer lower limit on Q by requiring the system to 
remain trapped in a resonance for a considerable fraction of its lifetime. Given that m2 = 0.067A/q 
and m2 ~ 10"^ Mq yr"\ the lifetime during which this system can sustain its high X-ray luminosity 
is estimated to be 7 million years, so a reasonable fraction of its lifetime to remain trapped in a 
resonance is at least 10^ yr. According to our model for (e/0. 009)^2/^2 > 4.0 x 10^ the system re- 
mains trapped in the resonance for more than 10^ yr (see Figure[8]), and as Figure |9] demonstrates, 
the mass transfer rate remains within roughly 10% of its nominal value m2 ~ 10"^ Mq yr"'. In 
this case, the eccentricity of the inner binary will never damp down to a fixed point because it is 
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indirectly driven up by semimajor axis expansion due to mass loss on a timescale of order 10"^ yr, 
which is at least an order of magnitude shorter then the timescale for eccentricity damping due to 
tidal dissipation. As expected from the evolutionary scenario the intrinsic period derivative is pos- 
itive, but because of the effect of tidal dissipation it is smaller than that due to Roche lobe overflow 
alone (see Figure [9l); the nominal value for the peri od derivative due to Ro che lobe overflow alone 
for our system parameters is P /P ^ 1.3 x 10"^ yr"^ (|Rappaport et al.ll 19871) . 

Finally, we note that if the inner binary is in fact expanding, the eccentricity will tend to in- 
crease as well. If the eccentricity is large enough, then Roche lobe overflow will occur through 
both the inner and outer Lagrange points, in contradiction with the low observed x-ray absorp- 
tion. Figure [8] shows that the eccentricity, while increasing with time, rem ains smaller than .07, 
consistent with the lack of mass loss through the outer (L2) Lagrange point (|Regos et al.ll2005() . 



5.1. The nature of the third body 



If the outer star is a white dwar f or a main sequence star , its mass is constrained to be < 0.5A/q 
by the lack of an optical detection dChou & Grindlayll200l|) . The lack of absorbing material along 
the line of sight to the X-ray source ind icates that the th ird star is not overflowing its Roche lobe. 
From the Roche lobe fitting formula of (|Eggletonlll983|) . 



R3<Rl 



0.49^2/3 



0.6^2/3 + in(l +^1/3) 



(29) 



where q is the mass ratio of th e third star to the to tal mass of the inner binary. This translates to 
7?3 < O.36i?0. From table 9 inl seattv et"aD (l2007h . this implies < 0.39Mq. We conclude that 
the only stars with mass > OAMq that will fit into the outer orbit is a white dwarf or neutron star 
(or black hole). This leaves open the possibility that the third star is a main sequence star with 
m < OAMq. 



According to llvanoval (l2008h . her Table 1, the fraction of hierarchical triples consisting of a 
neutron star primary, a white dwarf secondary, and a white dwarf tertiary formed via binary — binary 
encounters is about 1.4 x 10"^. The fraction of triples consisting of a neutron star-white dwarf bi- 
nary orbited by a O.4M0 (or lower) main sequence star is similar. The fraction of neutron star-white 
dwarf-neutron star systems is 2.1 x 10"^. If this is the primary channel for formation of triple star 
systems, the third star is likely to be either a white dwarf or a low mass (m < OAMq) main sequence 
star. 
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Fig. 8. — The eccentricity as a function of time (upper panel) and the argument of perias- 
tron as a function of time ( co versus t, lower panel) in the quadrupole approximation using 
(e/0. 009)^2/^2 = 4.5 X 10^. The system remains trapped in the resonance for more than lO^yr 
which is a considerable fraction of the system lifetime. The eccentricity stays under the limit of 
0.07, a constraint imposed by the absence of L2 mass loss. 
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line) as a function of time in the quadrupole approximation using (e/0. 009)^2/^2 = 4.5 x 10^. The 
system remains trapped in the resonance for more than 10^ yr, which is a reasonable fraction of the 
system lifetime. The mass transfer rate is within 10% of its nominal value mi w 10"^ M0 yr"^ 
P/P is lower than that due to Roche lobe overflow alone (dashed line), but still > 0. 
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6. DISCUSSION 

The origin of the 170 day luminosit y variations in 4 U 1820-30 was first attributed to the 
pr esence of a third body in the system by Grindlay J 19881): this possibi li ty was expanded upon 
by Ichou & Grindlavl J2OO1I) and more recently bv lzdziarski et all (l2007h . Izdziarski etaP J2OO7I) 
used a numerical model that calculates the time evolution of an isolated hierarchical triple of point 
masses, using secular perturbation theory up to octupole terms. Their model neglects the effects 
of tidal and rotational distortion of the white dwarf, tidal friction, mass transfer and gravitational 
radiation from the inner binary. Their calculations do include the GR periastron precession of the 
inner binary. 



Zdziarski et al.l (|2007|) find a configuration that reproduces the 171 day oscillations (assuming 



they are due to variations in e). They note that the GR precession rate is near 170 days, and then 
choose a rather low neutron star plus white dwarf mass of 1.29-I-O.O7M0. With this choice, the 
period of the GR precession is ~ 168 days. This period is very near, but slightly shorter than, the 
observed 171 day period. To arrive at the longer period, they choose the location and inclination 
of the third body so that the Kozai torque results in a retrograde precession. When added to the 
GR precession, this retrograde Kozai precession ensures the period of eccentricity oscillations will 
be longer than 168 days. They are driven to a much lower magnitude for the Kozai torque than 
employed in this paper; they use agut = 8.66a and iq = 40.96°. They start with tu = 0° and e = 10""^, 
ensuring that their solution circulates rather than librating. 

They note that the apparent near equality between the Kozai and GR precession rates is "a 
very remarkable coincidence", but go on to say that they do not have any explanation for this 
coincidence. 

We have argued that the origin of the 171 day period of the luminosity variation of LMBX 
4U 1820-30 arises from libration in the Kozai resonance. This trapping explains why the Kozai 
precession rate is comparable to the sum of the other precession rates in the problem. If ^2 is small 
enough, then the largest precession frequency in the absence of a third body is that given by general 
relativity. In that case, the Kozai and GR precession rates will sum to zero, i.e., the magnitude of 
the two precession rates will be equal. Hence if ^2 is small, then the e xpansion of the orbit o f the 
inner binary naturally explains the "remarkable coincidence" noted by lZdziarski et al.l (l2007b . We 
stress that, independent of the value of k2, the natural state of the system is likely to be libration 
rather than circulation. 

Trapping into libration is a consequence of mass-transfer driven orbital expansion in the inner 
binary. We have pointed out that the apparent negative period derivative, if it were intrinsic to 
the system, would not last for reasonable fraction of the system's lifetime. We find this to be an 
untenable situation. 



The observed negative period derivative of the inner binary allows us to constrain the tidal 
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dissipation factor Q yielding a very firm lower limit of (e/0. 009)^2/^2 > 3.15 x 10^. We argue, 
however, that (e/0. 009)^2/^2 has to be yet higher, to trap and maintain the system in libration 
around the stable Kozai fixed point. Our finding indicates that if 4U 1820-30 is indeed a triple 
system, the negative period derivative is not an intrinsic property of the system. However, as we 
showed in section [5] it does not arise from the accele ration of the gravitation fi eld of the globular 
cluster in which 4U 1820-30 resides, as suggested by Ivan der Klis et al.l (Il993b|) . 

In general, the eccentric orbit of a close binary system similar to 4U 1820-30 could lead to 
a time-dependent irradi ation of the secondary which could, in turn, give rise to a thermal tide 
(I Arras & Socratesll2010f) . A thermal tidal torque opposes the gravitational tidal torque, tending 
to force the secondary away from synchronous rotation and to enhance the orbital eccentricity. 
An asynchronous spin may cause large tidal heating rates, depositing heat in the interior of the 
secondary. In addition, the irradiation of the stellar surface by the neutron star (or by the accretion 
disk) will reduce the heat flux from the center of the white dwarf outward, so these irradiated white 
dwarfs will be hotter than passively cooling white dwarfs. Since they are hotter, they will have 
larger radii. The interplay between the two tidal torques would eventually set the equilibrium spin 
state. As long as this equilibrium state is not reached, the resulting bulge may oscillate, causing 
a periodic exchange of angular momentum between the orbit and the spin of the white dwarf. 
This might provide an alternate mechanism for producing the lumi nosity variations in 4U 1 820- 
30. Since this period is very stable, P^/P^ < 2.2 x 10""^ according to lChou & Grindlayl (1200 lb . we 
are currently looking into possibility of such an interplay between gravitational and thermal tidal 
torque as an explanation for 171 day period in 4U 1820-30. 

We anticipate that the resonance trapping mechanism we have described in this paper is 
generic in Roche lobe overflow binaries in triple systems. The exact nature of the librating or- 
bit will vary with the properties of the particular system. For example, for a binary with a larger 
semimajor axis, such that utb > tucR, resonant trapping will lead to |(u)fa,2a, | = logr- 



7. CONCLUSIONS 

This paper provides an estimate for a lower limit of the tidal dissipation parameter Q for a 
Helium white dwarf. It also elucidates the possible evolutionary history of 4U 1820-30, i.e., how 
the system arrived at a state where the secular dynamics are not dominated by the effects of the 
white dwarf's tidal bulge, despite the fact that the white dwarf is overflowing its Roche lobe in an 
orbit with a period of 685 s. 

We suggest that the system is trapped in the Kozai resonance. This resonance trapping is 
responsible for the observed 171 day period, which we interpret as the period of small oscillations 
around a stable fixed point in the Kozai resonance. If the system is not librating, one requires very 
fine tuning to get the 171 day period. 
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We provide lower limit on the tidal dissipation rate, as measured by the factor Q; (e / 0.009)^2/ ^2 > 
4x lO''. 

Further exploration of the long term (tidal and mass overflow-driven) evolution of this and 
similar short period ultra compact X-ray binaries is clearly warranted. Inclusion of the thermal 
tides into dynamics of these systems may introduce an alternative explanation for origin of long 
period modulation of the light curve. For the particular case of 4U 1820-30, better modelling of 
the gravitational potential in the host globular cluster, NGC 6624, would allow for an upper limit 
on Q. We are pursuing both lines of investigation. 

The authors are grateful to Natasha Ivanova, Cole Miller, Doug Hamilton, Andrew Gumming 
and Phil Arras for helpful discussions. This research has made use of the SIMBAD database, 
operated at CDS, Strasbourg, France, and of NASA's Astrophysics Data System. The authors are 
supported in part by the Canada Research Chair program and by NSERC of Canada. 



APPENDIX A 
EQUATIONS OF MOTION 

The equations of motion we employ model the Kozai interaction, the dynamical effects of the 
tidal bulge of the He white dwarf, GR periastron precession, the rotational bulge of the He white 
dwarf, conservative mass transfer driven by the emission of gravitational radiation, and tidal dis- 
sipation. We do not consider tides raised on the neutron star primary. Detailed derivation of the 
e quations representing K oza i cycles with tidal friction and GR per iastron precession can be found 



in 



Eggleton et al.l (|1998|) and lEggleton & Kiseleva-EggletonI (1200 Stellar masses are denoted by 



mi (the mass of the neutron star primary), m2 (the mass of the white dwarf secondary) M = m\+m2 
(the inner binary mass), nii, (the mass of the outer companion), and the reduced mass of the inner 
binary ^ = mim2/(mi+m2). The mean motion of the inner binary is n = 2Tc/P= [GM/a^y^^. The 
inner binary orbital elements are: semimajor axis a, eccentricity e, mutual inclination between the 
inner binary and the outer binary orbit i, the argument of periastron to, the longitude of ascending 
node ri. is the tidal Love number, Q is the tidal dissipation factor, and is the radius of the 
white dwarf. The orbital parameters of the outer binary are denoted Uuut and Cout- G is Newtons 
constant and c is the speed of light. 

Changes in the semimajor axis of the inner binary a are caused by tidal dissipation and mass 
transfer: 

d = drD + aMT, (30) 
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For the zero eccentricity case, equation [34| is derived in detail in lRappaport et al.l (I1987|) . where in- 
stead of the dependancy on periastron period Pperiastmn they consider dependancy on binary period. 

The tidal friction time scale is 
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(35) 



The eccentricity of the inner binary e is affected by the Kozai torque and by tidal dissipation: 



where 
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The mutual inclination between the inner and the outer binary orbit, i is affected by Kozai 
torques, by the rotational bulge, and by tidal dissipation: 



I - iKozai + lRB + lPD, 



where 
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(1-^2)5 



(40) 

(41) 
(42) 
(43) 
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Besides the negative precession rate of the argument of periastron due to Kozai cycles, the 
total precession rate of the argument of periastron has additional positive contributions from the 
tidal bulge, GR, the rotational bulge, and the tidal dissipation: 



■ ^ Kozai + '^TB + '^GR + ^RB + 1^70, 



where 



Gm3 



UJtb — 



^GR - 



(^RB - 



X 



15(GM)5 8+12^2+^4^1 



16a 2 
3(GM)i 

a^c\l-e^) 



2(l-e^) + 5 sin^ uj(e^ - sin^ i) 



M-2 k2Rl 



4G2a2(l-e2)2 mz 
(lill-^ll-il^) +2f2/iCOt« sine; + 0^ cos 



^TD — ~ 



QeCOs(oj)coU \ + \e^ + \e^ ri^sin(a;)cotj \-^\e^ + \e^ 



(l-e2)5 



■ + ■ 



(l-e2)5 



(44) 

(45) 

(46) 
(47) 
(48) 
(49) 
(50) 



The precession of the longitude of ascending node is caused by Kozai cycles, rotational bulge and 
tidal dissipation: 

^in = ^Kozai + ^RB + ^TD, 



where 



^Kozai 
^RB 



COSi 



3 + 12e2-15e2cos2a; 



mi ^2^2 



- n„ COS a; - n„ sin w 



l/jina^ sin? {\-e^Y 
fie cos a; 1 + + 1^"^ Qq sin a; 1 + + fe"^ 
Insimtp (l-e^)^ Inmntp (l-e^)^ 



(51) 

(52) 
(53) 
(54) 
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APPENDIX B 
ADIABATIC INVARIANCE OF THE ACTION 



Time-dependent Hamiltonians, even those with just one degree of freedom, can be difficult to solve. 
However, for Hamiltonians where the time dependance is sufficiently slow, the problem is easier to 
tackle due to the existence of variables that are almost constant. The approximate constants are the 
action variables of the Hamiltonian, when the slow time dependance is neglected. Suppose that the 
time dependance enters through a time dependent parameter K{t). If the parameter k varies very 
slowly with time, treating n as time-independent parameter allows us to find action-angle variables 
following the standard prescription. These action-angle variables are function of time through K{t), 
which leads to the action no longer being a constant of motion. However, when k varies slowly 
with time, the action is nearly constant. Such an action is known as an adiabatic invariant. 

As described in section |3l capture in the resonance is a natural consequence of semimajor 
axis expansion driven by mass transfer from the white dwarf. The Hamiltonian of our system 
(see equation IH) is a function of the semimajor axis, which is a parameter of H, playing the role 
of K(t). In our case the semimajor axis is not the only parameter varying with time; the masses 
of the inner binary vary with time as well. Here we show, both analytically and via numerical 
integration, that the change in the eccentricity is coupled to the change in the semimajor axis. 
When the semimajor axis expands (respectively, contracts) the eccentricity of the stable fixed point 
increases (decreases). We also demonstrate that the timescale for the change in the eccentricity is 
a factor of > 150 shorter than the timescale for the semimajor axis. 

To find the action we expand our Hamiltonian (equation ID) around the fixed point: 



Since we are expanding around the resonance, all terms cx A ^ vanish. After some algebra we find: 



which is similar to the Hamiltonian of the harmonic oscillator. Written more compactly (and 




(55) 
(56) 
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implicitly defining a(t), j3{t) and C{t)): 



(58) 



We solve for A Q and evaluate the integral 



71" ./n V a a 



/o V « 
where /\u,„ax = {2(Ho- C{t))/(5) ^ . We find: 

Ho-C{t) 



(59) 



(a/3)^ 



7: 

Plugging in the corresponding terms from equation (1571) yields: 

C(a,mi,m2) Pi(e/,a,mi,m2,«/) 



(60) 



J: 



e/smif p2 



(61) 



P^(ef,a,mi,m2,if) 



where 



— - 1 - 1 2( 1 - e^.) cos^ If + 9( 1 - e^) + 1 5 cos^ + - ( 1 + ) + 4/^2 -(2 + 1 Se^.) 
A A2 A 



D 3 , 
+^2^(1 + 2^) 



and 



^2 = 30 



B 



D 



1 8 - 24cos' + 2-(l + -e)) + 120yt2-(13 + 84e}) + 12/t2-(l + -4). 



Since the action J is an adiabatic invariant, we have: 

dJ dJ 



dJ _dJ . dJ . 

dt de da \ dm2 dnii 



m.2 = 0. 



The partial derivatives are: 



dJ 


J 


de 




dJ 


- -i 


da 


a \ 


dJ 


J 


dm2 


m2 


dJ 


J 


dm I 


nil 



Sf V Pi de 2P2 de J ' ef 



a dPi _l a dP2 
Fi~da 2^aa" 
m2 dPi 1 m2 dP2 



a 



1 + 



Pi dm2 2 P2 dm2 
nil dPi 1 nil dP2 



= C — 

' ni2 

Pi dnii 2 P2 dnii I ' nii 



(62) 
(63) 

(64) 
(65) 

(66) 
(67) 
(68) 
(69) 
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(72) 

a 



Plugging these partial derivatives back into equation (1651) yields: 

= C- + a- + C,n, -Cn,A —. (70) 

Cf a \ nil J m2 

The inner binary orbit is eccentric, which makes the mass transfer rate proportional to perias- 
tron distance = a(l -e). Hence the m2 term can be decoupled into two terms, one proportional 
to e and the other proportional to d: 

^ = ^Ji = -l^ + l-^ (71) 
m2 rp la 2\-e 

Combining equations (TTOl) and dTTI) and solving for e leads to: 

1 ( C -^C ^ - C 

e 2\^'"2 mi ^'"i J ^" 

Cf n I 3 g/ In _Jn2r 

Plugging in the numerical values: 

e d 

— ^ 150-. (73) 

ef a 

Defining the time scales for the eccentricity and the semimajor axis to decay or increase (depending 
on the value of Q) as Te = ef/e and = a/d, the timescales in equation (1731) are related by: 

Te ~ 6.7 X 10"V«. (74) 

To demonstrate that the eccentricity evolution is indeed a consequence of the action being an 
adiabatic invariant, we follow the evolution of the orbit around the fixed point Cf = 0.0155 and 
ujf = 90°. Figure [To! shows e as a function of time in a case where the semimajor axis is increasing, 
meaning that tidal dissipation is sufficiently weak so that the evolution of the semimajor axis is 
dominated by mass transfer (2 = 8 x 10^). The solid line presents e predicted by equation ([731) . 
For t > 10^ yr the numerical integration gives e ~ 10"^ yr"' corresponding to a timescale 150 times 
shorter than the timescale for the semimajor axis. 

Despite the fact that the semimajor axis is expanding, the numerical integration shows a tran- 
sient phase (roughly the first 2000 years) where de/dt < 0, and a longer phase (~ 10^ yr) where 
e is larger than predicted by equation (1731) . There are contributions to the eccentricity evolution 
which we have ignored in our analytic treatment; for example the spin of the white dwarf is not 
locked during the evolution of the system. These un-modelled contributions are the source of the 
transient behaviour. 

To support this statement, we illustrate the eccentricity evolution in various cases where we 
turn off different dynamical effects in Figure \TT\ The solid line presents a result from the numer- 
ical integration that includes all dynamical effects in our model, while the dotted line is the same 
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105 2x105 3x105 4x105 5x105 

t[yr] 

Fig. 10. — e as a function of time for the case where the semimajor axis is expanding, 2 = 8 x 10^. 
We start integration exactly at the fixed point, where initial eccentricity is e/ o = 0.01555. All other 
parameters are as listed in table 1 . The solid line comes from the analytic estimate where the action 
/ is considered to be an adiabatic invariant. The dashed line is a result of numerical integration. As 
expected from the action being adiabatic invariant, e is positive. The difference in the magnitude 
of e within first 2 x 10^ yr is a result of our simplified analytic calculation that does not include spin 
dynamics. Since our analytic estimate is valid for small eccentricities, here we stop the integration 
when e > 0.1 
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integration with the ejo term set to 0; the result shows that direct tidal dissipation on the eccen- 
tricity (equation [381) is not dynamically significant. The dashed line presents the case where djD is 
set to (see equation [3T|); The result shows that ajD has a significant influence on the eccentricity 
evolution. The dash-dotted line shows the eccentricity when tidal dissipation factor Q is set to 
infinity, but only in the differential equations that govern spin evolution. The long dash-dotted line 
shows the eccentricity evolution in the case where Q is set to infinity in the equations that govern 
the evolution of the spins together with aVo being set to 0. The latter three cases demonstrate that 
the eccentricity starts increasing immediately with the semimajor axis expansion, which is exactly 
the behaviour predicted by the analytic analysis. After 5.5 x 10^ yr the eccentricity becomes > 0. 1 
and since our analytic estimate is valid only for small eccentricities we stop the integration here. 

The cause of the transient behaviour is the fact that the spin of the white dwarf is not, contrary 
to our choice of initial conditions, tidally locked during the evolution of the system. Whether the 
spin settles down in some Cassini state or other stable configuration later during the evolution of 
the system is a possibility open to further investigation. 
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2x10'' 4x10'' 6x10'' 8x10" 10^ 

t[yr] 

Fig. 11. — The eccentricity of the fixed point of the inner binary as a function of time. Initial 
conditions are the same as in Figure \T0\ The solid line is a result of the numerical integration 
including all dynamical effects, while the dotted line is a result of the same integration but with 
eVo = 0; these two results demonstrate that direct tidal dissipation on the eccentricity does not 
significantly affect the evolution of eccentricity. The dashed line is the result of integration where 
ciTD = 0; this term has a more significant effect on the evolution of the eccentricity. The dash-dotted 
line presents the case where we set Q = oo, but only in the equations that govern the spin evolution. 
The long dash-dotted line shows the eccentricity evolution when Q = oo for the spins and dro = 0; 
in this case we have the fastest increase in the eccentricity. 
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TABLE 1 . System parameters 


Symbol 


Definition 


Value 


Reference 


nil 


Neutron star (primary) mass 


IAMq 




m2 


White dwarf (secondary) mass 


0.067Mq 


RaDDaDort et al. ("1987) 


m-i 


Third companion mass 


0.55Mq 




ai 


Inner binary semimajor axis 


1.32X lO^^cm 


Stella etal. (1987) 


^out 


Outer binary semimajor axis 


S.Oai 




^in.O 


Inner binary initial eccentricity 


0.009 




Sout.O 


Outer binary eccentricity 


10-4 




^init 


Initial mutual inclination 


44.715^ 






Initial argumet of periastron 


90" 






Longitude of ascending node 







Ri 


White dwarf radius 


2.2 X 10"^ cm 




ki 


Tidal Love number 


0.01 


Arras (private communication) 


Q 


Tidal dissipation factor 


5x 10^ 
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